Measurement of Effective Temperatures in an Aging Colloidal Glass 
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We study the thermal fluctuations of an optically confined probe particle, suspended in an aging 
colloidal suspension, as the suspension transforms from a viscous liquid into an elastic glass. The 
micron-sized bead forms a harmonic oscillator. By monitoring the equal-time fluctuations of the 
tracer, at two different laser powers, we determine the temperature of the oscillator, T . In the 
ergodic liquid the temperatures of the oscillator and its environment are equal while, in contrast, in 
a nonequilibrium glassy phase we find that T substantially exceeds the bath temperature. 



Understanding the slow dynamics of glasses is one of 
the most fascinating yet difficult challenges in statisti- 
cal physics. One question which has attracted consider- 
able interest is whether the dynamical fluctuations of a 
glass can be characterized by a non-equilibrium temper- 
ature. A rigorous thermodynamic temperature is strictly 
impossible to define for an aging glass, which remains 
far from thermal equilibrium even on long-time scales. 
Nevertheless various groups [J, 0] have proposed "ef- 
fective" temperatures with many of the properties ex- 
pected for an equilibrium temperature. These ideas lead 
to the surprising prediction of two distinct temperatures 
in an aging glass [lj . The fast rattling modes of parti- 
cles inside the cage constituted by neighbors thermalizc 
rapidly to the temperature of the environment, Tbath, 
while the much slower structural rearrangement of these 
cages are supposed to be characterized by a second tem- 
perature, T g, which mean-field models predict should 
exceed Tbath- To date most of the support for this strik- 
ing two-temperature picture has emerged from simula- 
tion results [2| on idealized glasses. Experiments have 
so far produced conflicting results. Studies of colloidal 
glasses have reported that T e g increases remains un- 
changed Q , or even decreases (HI with age in contrast 
to measurements on structural |(| and spin glasses Q 
which have revealed effective temperature warmer than 
the bath temperature. 

In this Letter, we report the temperature of a 
micrometer-sized sphere immersed in an aging colloidal 
suspension, as the suspension transforms from a fluid to 
a glass. The particle, captured in an optical trap, con- 
stitutes a microscopic harmonic oscillator whose fluctu- 
ations probe the nonequilibrium dynamics of the aging 
glass. We measure the equal-time fluctuations of this 
local oscillator and show that, in an ergodic phase, the 
temperature of the oscillator T Q equals the environment 
temperature Tbath of the system, as required by equilib- 
rium statistical mechanics. Significantly, when we repeat 
the measurements in an aging glass we find a higher tem- 
perature and T Q > Tbath- 

All experiments to date on colloidal glasses 0, 0, Q 
have relied on active, driven measurements. We conduct 
our experiments instead in a quasi-static limit which has 
several advantages. First, the experiments are simpler 



because there is no need to characterize the complete 
time-dependent response of the system. Second, the lack 
of an external driving force ensures that we always op- 
erate within the linear response regime. We determine 
the temperature of the laser-trapped particle by recall- 
ing that the fluctuations in the coordinate, 5x = x — (x) , 
of an oscillator (with spring constant /ct), in contact with 
an equilibrium system at a temperature T, are Gaussian, 
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Here fcs — 1 and the static susceptibility \ ~ the shift 
in the mean position (x) produced by a small constant 
external field - is simply l/hr- Operationally, we define 
the temperature of the oscillator, T Q , by the ratio between 
the one-dimensional mean-squared displacement (MSD) 
of the probe and the static susceptibility, T Q = (fa 2 ) j\. 
Equilibrium statistical mechanics guarantees T Q = T, ir- 
respective of the nature of the coupling between oscillator 
and system. Out-of-cquilibrium however this equivalence 
no longer holds. Berthier and Barrat [1] have argued that 
a generalized cquipartition principle holds in nonequilib- 
rium glassy materials, implying that an oscillator will 
record the (slow) effective temperature of a glass. We 
postpone the interpretation of T Q until later and use it, 
for now, simply as a experimentally accessible measure 
of the local temperature of a material. 

Brownian fluctuations change considerably at the glass 
transition simply because of elasticity. A fluctuating 
optically-trapped sphere in an elastic medium, such as a 
glass, is subject to two harmonic forces; an optical force 
Topt, resulting from the external laser field, and an addi- 
tional elastic force F c , caused by matrix deformation. In 
an infinite continuum, with shear modulus G and Poisson 
ratio v, the elastic force is a linear function of position, 
F c = —k c (x — (x)). The spring constant is k c — 6irGR 
for a sphere of radius R, when v — 1/2 The optical 
forces are T opt = — k[x — (x)), with k the optical trap 
stiffness, so the total spring constant of the trapped par- 
ticle is &t = k c + fc. To separate the elastic contribution 
from the optical contributions to fcx we modulate the in- 
tensity of the laser beam, rapidly switching it between 
two levels, which we identify below by the subscripts 1 
and 2 respectively. Since fc is independent of laser power, 
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consecutive measurements of (Sx 2 } at two different laser 
powers provides estimates of both k e and T Q , 

k 2 (Sx 2 ) 2 -k 1 (Sx 2 ) 1 

{8x 2 ) l -{8x 2 ) 2 
{k 2 -k l ){5x 2 ) 1 {5x 2 ) 2 
°~ (6x 2 )^(8x 2 ) 2 ' () 

Here fej (i = 1, 2) is the optical force constant at each 
laser power. 

The glass studied was a transparent aqueous disper- 
sion of charged colloidal disks of radius a w 15 nm and 
1 nm thick (Laponitc RD). The suspension (2.4% mass 
fraction) was filtered through a 0.45 /im filter to obtain 
a reproducible initial liquid at t = 0. Immediately af- 
ter filtration, a small amount of a dilute suspension of 
silica spheres (radius R = 0.55 ± 0.03 /im) was added. 
A particle was captured in a tightly-focused laser beam 
(A = 1064 nm) and its lateral position, x(t), measured 
with a quadrant photodetector. The strength of the opti- 
cal trap was cycled every 52 s between an initial (optical) 
trap strength, k\ = 4.4 pN//im, and a final trap strength, 
fc 2 = 11.0 pN//im. All experiments were performed at 22 
± 1 °C. We confirmed that laser heating was negligible 
by following the temperature of a trapped silica bead 
in deionized water. Recording continuously for over an 
hour at ambient conditions gave a mean temperature of 
T = 300 ± 8 K 0, with no upward drift in T . Re- 
peating the measurements on different particles showed 
that the small variation in bead size led to a systematic 
error in T Q of about ± 40 K. 

After mixing the suspension gradually thickens with 
time and becomes jammed in a glassy state [llj]. To 
identify the point of dynamical arrest we measure the 
autocorrelation function of the Brownian fluctuations of 
the probe sphere, C x (r) = (Sx(t)Sx(O)) / (<5x 2 ). The po- 
sitional fluctuations are well represented by the stretched 
exponential relaxation, C x (r) = exp [^(t/tq)' 3 ] , familiar 
from generic glassy systems, where r Q is a terminal re- 
laxation time and (3 is a stretching exponent. The inset 
to Fig. [1] reveals that the fitted exponent j3 varies con- 
siderably with age. For t < t c (« 220 ± 15 min) the 
correlation function is stretched (/3 < 1) but with in- 
creasing age there is a crossover from a stretched- to a 
near-exponential form. To decouple changes in the re- 
laxation time from the changes in the exponent (3 we 
evaluate the 'mean' relaxation time, r m = (r Q //3)r(/3 _1 ) 
where V is the gamma function. The relaxation dynam- 
ics clearly exhibit two distinct regimes. There is first a 
fast regime where r m grows approximately exponentially 
with sample age, followed by a second regime of much 
slower growth where, T m ~ t. Dynamical arrest occurs 
at a time t c which we estimate from the intersection of 
the corresponding fits in Fig. [T] as 220 ± 15 min. For 
t < t c , the system behaves as if it were an ergodic liquid, 
albeit one with a relaxation time which slows down pro- 
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FIG. 1: Variation of the mean relaxation time r m with time 
after preparation for hi — 11.0 pN/ fim. The correlation func- 
tion C x (t) was fit to a stretched/compressed-exponential re- 
laxation. The continuous lines are fits to the exponential and 
linear aging regimes respectively. t c identifies the crossover 
from exponential to the full aging regime and indicates the 
approximate position of dynamical arrest. Similar plots are 
found if the relaxation time is chosen as r Q or defined as the 
time at which C x (t) has decayed by a factor of 1/e. Inset: 
The age dependence of the fitted exponent /3. 

gressively as t increases. At t « t c the suspension arrests 
and begins to age. 

To explore the effect of the glass transition on the 
temperature recorded by the Brownian probe, we first 
examine the probability distribution of the fluctuations, 
P(5x). Fig. reveals that although the particle is 
highly constrained by the glassy matrix, the trajectories 
exhibit fluctuations which remain Gaussian. This implies 
that the Einstein expression for the fluctuations (Eq. [T|) 
can be generalized by replacing the bath temperature 
with an effective temperature. To measure the effective 
temperature we determine the MSD of the trapped par- 
ticle, recognizing that the dramatic dynamical slowing at 
the glass transition, evident in Fig [TJ requires care if the 
asymptotic limit is to be evaluated. To check for con- 
vergence, we split each recorded 52 s particle trajectory 
into a number of shorter equal-time lengths of duration 
ts and evaluate the MSD of each, before averaging the 
values together. We choose is = 3.3 s in order to have 
a good signal to noise ratio whilst retaining many long- 
lived fluctuations. The elastic modulus G of a (hard- 
sphere) repulsive glass scales as G ~ fc^T/a 3 , where a 
is the radius of the particles that form the glass, so the 
thermal fluctuations of a sphere of radius R embedded in 
a glass are restricted to a cage with a characteristic size 
of l g ~ (a 3 /^) 1 / 2 0, Taking the effective radius as 
a ~ 45 nm [13[ yields I 2 « 200 nm 2 for our system. Con- 
sistent with this picture the asymptotic MSD exhibits 
two distinct regimes. Fig. [D(b) demonstrates that for 
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FIG. 2: (color online), (a) Probability distributions of the 
fluctuations of probe particles for different times after prepa- 
ration, t = 200 min (circles), 250 min (triangles), and 301 min 
(squares). The dashed lines are Gaussian fits to the data, (b) 
Asymptotic mean-square displacement (MSD) of probe, con- 
fined in an optical trap of strength k\ = 4.4 pN//im, as a 
function of sample age t. The statistical errors are of order 
±11 nm 2 . The arrow indicates the square of the localization 
length in the glass, lg. The solid line is an exponential fit to 
the MSD within the glass. The inset depicts a linear plot of 
the MSD (in nm 2 ) as function of age, in weak (fei: circles) 
and strong {k%: squares) laser traps, k% = 11.0 pN//im. Note 
the MSDs are identical, within error, for t > 250 min where 

<^ 2 > < z 2 . 



t < 100 min, where (dx 2 } > l 2 , the degree of localiza- 
tion does not change with age although it depends sensi- 
tively on the laser intensity, as expected for confinement 
by an external field. With increasing age, by contrast, 
the asymptotic MSD shows a progressively weaker depen- 
dence on intensity. Indeed when (5x 2 ^ ~ I 2 the thermal 
fluctuations are almost totally independent of the laser 
field, consistent with permanent caging of the probe par- 
ticle by the glass. The localization length in the glass 
shrinks exponentially as the glass get older, until at very 
long times (i ~ 350 min) all measurable motion ceases. 

The change in localization with time, evident in 
Fig. HJb), allows us to monitor the temperature of the 
slow structural modes of the glass. Figure [3] shows the 
age dependence of the oscillator temperature T Q and the 
shear modulus G, calculated from Eq. [2] and the data 
of Fig. Elb) (assuming v = 0.5). The data reveals a 
marked increase in T approaching the liquid-glass tran- 
sition. For t < t c the temperature of the Brownian probe 
is approximately constant at 247± 20 K [llj , which given 
that the systematic uncertainty in T is at least 40K, is 
in reasonable agreement with the environment temper- 
ature, Tbath = 295 K. By contrast, in the glass we find 
a systematic increase in the oscillator temperature with 
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FIG. 3: (color online) (a) Effective temperature as a func- 
tion of sample age. The circles and squares depict the re- 
sults of measurements using similarly-sized particles in two 
physically-distinct laponite samples. The scatter illustrates 
the typical level of reproducibility achieved. Errors bars, 
which are estimated from the variability in (<5a; 2 ) for t < 100 
min where the MSD is constant, are only shown on a few 
representative points for clarity, (b) Age dependence of the 
single particle elastic spring constant k e and the equivalent 
elastic modulus G = k e /6nR. 



increasing age, with the ratio T /Tb a th ~ 2 at the longest 
accessible age. Since the glass is non-ergodic an ensem- 
ble average quantity such as the effective temperature 
strictly requires an average over many probe particles. 
To verify the reproducibility of our results we repeated 
measurements using similarly-sized particles in different 
physical samples with the same nominal laponite con- 
centration and found no significant difference in the age 
dependence of T, /Tbath. An example of the data spread 
for data sets recorded from two different samples is shown 
in Fig. [3Ja). We attribute this homogeneity to the fact 
that the probe is significantly larger than the laponite 
particles. 

The growth evident in T Q as the glass ages is, at first 
sight, very surprising. Intuitively, one expects a glass 
to cool as it ages and T e ff to fall with increasing t. To 
account for the temperature increase, we analyze a sim- 
ple stochastic model of our experiments, ignoring the de- 
tailed microhydrodynamic coupling between the fluctu- 
ations of the embedded probe and the sea of platelets. 
The model, introduced in [l4[ , considers the diffusion of 
a harmonically-bound particle coupled to two thermal 
baths. The two baths are kept at different temperatures 
to mimic the effects of the fast and slow modes of the 
glass on probe diffusion. The fast bath is held at a tern- 
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FIG. 4: (color online) Elasticity dependence of the scaled 
oscillator temperature for T e s = 15Tb a th- The solid lines 
are obtained from the analytic solution of Eq. [3] for relative 
timescales = 10 (red), 100 (green), and 1000 (blue). 

Inset: Experimental dependence of To/Tbath on the elasticity 
ratio, fce/fe (circles). The continuous line depict the theoret- 
ical predictions for (, s /£f = 100. 



an aging glass. Using a generalized equipartition princi- 
ple we determine the temperature of the oscillator formed 
by the trapped probe. In a glass we find an oscillator tem- 
perature which is substantially higher than the temper- 
ature of the environment; while in the fluid the temper- 
atures of the oscillator and environment are essentially 
equal. We propose a simple theoretical model for these 
noncquilibrium experiments which reveals that, at high 
elasticities, the probe thermalizcs to the effective temper- 
ature of the slow modes of the glass. Our findings agree 
broadly with the results obtained recently by Abou et al. 
and Strachan et al. jjjj], who used different techniques 
to measure T e g . Taken together this agreement provides 
strong experimental confirmation for the existence of an 
elevated effective temperature in a glass. Future work 
will look at using comparably-sized probe and glass par- 
ticles to study the role of dynamical heterogeneities on 
the effective temperature. 



perature Tbath and exerts an instantaneous friction of 
£f£(t) on the particle while the second, slower bath is 
maintained at a (higher) temperature T e g and is associ- 
ated with the memory function, T(t). In the overdamped 
limit, where acceleration can be neglected, Brownian dif- 
fusion is governed by the generalized Langevin equation, 

kx{t) = - [ dTT(t-T)x(T)-Cfx(t)+0(t)+l(t). (3) 
Jo 

The random forces due to the fast bath are Gaussian with 
( 7 (t)) = and (7(4)7(7-)) = 2£ f T hath 5(t-T), whereas the 
slow bath is characterized by (0(t)0(T)) = T c sY(t — t). 
Approximating the glass by a Maxwell fluid, with zero 
shear viscosity rj s and shear modulus G, the retarded fric- 
tion on a sphere of radius R is T(t) = k e exp(— k e t/^ s ), 
with k e = 6nGR and £ s = 6irr] s R 0]. Following 
Ref. fl4j , an analytic expression may be written down for 
the reduced oscillator temperature T = T /Tbath- Fig.Q] 
shows the dependence of the oscillator temperature on 
the elasticity ratio k e /k. Clearly while the oscillator tem- 
perature always lies between the temperatures of the fast 
and slow baths, the exact value is sharply dependent on 
k c /k. Only for k 3> k does T approach the temper- 
ature of the slow bath. Plotting the experimental data 
for T as a function of k e /k2 reveals a very similar trend 
(shown in the inset of Fig. 0|. Consequently we attribute 
the increase in the measured temperature with age to 
the growing strength of the mechanical coupling between 
particle and glass. 

In summary, we have measured the quasi-static fluctu- 
ations of an optically-bound probe particle immersed in 
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